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Third order Coulomb correction to ti threshold cross section * 

Y. Kiyo t 

Institut fiir Theoretische Physik E, RWTH Aachen, 
D - 52056 Aachen, Germany 

We report on our result of third order Coulomb correction to the cross section a{e^e~ — > tt) near threshold. 
Analytic expression for the Coulomb energy and wave function at the origin are obtained. We discuss the 
significance of the Coulomb correction to the threshold cross section and heavy quarkonium phenomenology. 



1. Threshold cross section 

The ti threshold cross section has the following 
schematic form in the conventional perturbative 
expansion of as, 

(T(e+e- ^t<)-aBorn[l + — + (—)' + •••], (1) 

V \ V y 

where UBorn is a Born cross section, and v = 
\/l " ^rrif/s is speed of produced t, i. The combi- 
nation of (as/v)" appears to all order in the per- 
turbation theory, known as Coulomb singularity. 
The perturbative expansion in is not applica- 
ble to the threshold cross section because v is of 
order as near the tt threshold ^/s ~ 2mt, and 
the Coulomb singularity (as/v)"' dominates the 
cross section. To obtain meaningful cross section 
the Coulomb singularity has to be summed to all 
order in as ■ The physical origin of the Coulomb 
singularity is instantaneous gluon exchange be- 
tween t and t which has a small spatial momen- 
tum \q\ ^ V 0. This is called potential gluon 
because of its propagator 



~ iTTCpas 

Vc{q) = —2 — 



(2) 



where Cf — 4/3, which is the Coulomb potential 
Vc = —Cpoislr in coordinate space. 
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This understanding leads to a quantum me- 
chanical description which is equivalent to resum- 
mation of the Coulomb singularity. People had 
been looking for an effective field theory (EFT) 
which makes the resummation systematic. Fi- 
nally a non-relativistic version of QCD was de- 
rived, called pNRQCD/vNRQCD EE). The EFT 
makes the resummation systematic based on non- 
relativistic power counting z; ~ <Si 1 in La- 
grangian level. Nowadays we understand the re- 
summation using the EFT, and higher order cor- 
rections are taken into account in the EFT frame- 
work systematically. The NNLO calculation for 
the total cross section was completed by sev- 
eral groups and now we are working on the 
NNNLO corrections using the EFT. 

The EFT classifies the corrections into three 
classes: 

• Hard corrections included in Wilson coeffi- 
cients of (composite) operators. 

• Potential corrections, which are non-local (in 
space) 4-Fermi operators but local in time. 

• Dynamical gluon corrections called ultra-soft 
gluon in the EFT. 

The hard corrections are related to re- 
normalization of the operators in the EFT. The 
ultrasoft correction appears at NNNLO calcu- 
lation, and the corresponding energy level cor- 
rection is know by Kniehl-Penin [5]. However 
complete NNNLO correction to the total cross 
section is not known, yet. We discuss a part of 
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the potential corrections in this report, which has 
the following form in the EFT Lagrangian, 

S£{x) = J d^r [V-^ V] (a; + r)V{r) [x^x] (^) , (3) 

where , X is a creation operator of heavy quark 
and anti-quark, respectively. We will report on 
the result of the Coulomb correction. 

We parameterize the momentum space 
Coulomb potential as follows 

V{q) = Vc{q) [1 + Vi{q) + ^q) + Vsiq)] , (4) 

where Vn{q) is n-th order correction to the pure 
Coulomb potential Vc , induced by loop diagrams 
when the EFT is derived from QCD. Explicit 
form of the potentials were derived at NNNLO 
(except 03) in ref. [HI, and the Coulomb potential 
reads 

^2 = (^)' [«2 + (2ai/3o + A) Ig + Pi ll ] 

^3 = (^)' h + 8^2c3 ;^(^,) + (302^0 + 2ai/3i 

+P2) Iq + (3ai/3o' + Ig + ] : (5) 

where Iq — ln(^^/q^) and Zq(i^) ~ ln(i^^/q^), f3i 
are the coefficients of QCD /3-function. The scale 
fi is QCD renormalization scale, and the v in 
V3 is a scale introduced to separate the ultra- 
soft and potential modes in the EFT. Physical 
quantities are scale independent if all the correc- 
tions at given order are taken into account (see 
for instance 7 ) . Now our task is to calculate the 
threshold cross section using the EFT Lagrangian 
eq.® with the Coulomb potential cqs.Q), ((S)). 

2. Green function method 

We use the optical theorem to calculate the 
threshold total cross section which tells us that 
the cross section can be obtained by taking the 
imaginary part of the correlation function of pro- 
duction currents 

a{e+e- ^ tt) = l^ImG(^ + iTt), 
mi 

G=^J rf'xe*«^(vac|[^t^](^)[^t^](o)|vac),(6) 



where Nc = 3 and mt is the quark pole mass. 
Using the EFT one can show that the matrix el- 
ement G{E) can be expressed by the quantum 
mechanical Green function 

G{E) = (0|G(i?)|0) = (0| — |0) , (7) 

where |0) denotes a quantum mechanical posi- 
tion eigenstate at the origin r = 0. At this 
stage we perform quantum mechanical perturba- 
tion theory by expanding the denominator of the 
Green function with respect to the higher order 
Coulomb potentials. The third order corrections 
read 

S^G = - GoSV-iGo + 2Gq5ViGq5V2Go 

- GQ5ViGo5ViGo5ViGa ■ (8) 

Here Go — {p^ /mt — E)~^ is the zeroth order 
Green function, 5Vn = [V„Vc'](r) is the n-th or- 
der Coulomb potential. We calculated the ex- 
panded Green function semi-analytically and ob- 
tained double sum representations, which were 
evaluated numerically p. 

The Green function G{E) has a single pole at 
the bound-state energy level E — En, 

G(E) = h non-smgular , (9) 

En — E " it 

where |V'n(0)P and En is the bound-state wave 
function squared at the origin and energy level, 
respectively, which has series expansion in as 

3 

|^„(o)P = l^i°no)P[i + EO^^^]' 

n 

where En\\'4'n(0)\'^ is the zeroth order result 
for the energy and wave function. By perform- 
ing matching between expanded Green function 
eq. © and the pole structure of exact Green func- 
tion eq.©, we obtained analytical result I] for 
Ci and fi for the S-wave bound state at NNNLO. 
The expression is too lengthy to show here. In the 
next section we discuss phenomenological signifi- 
cance of the NNNLO Coulomb corrections using 
the obtained result. 
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3. Numerical analysis and Conclusion 

Our formalism up to now is not specific to the 
top quark, actually it is applicable to bottom 
quark system replacing mt and nf = 5 by 
and nj = 4 (n/ exists in the as and /3-function 
in the Coulomb potential.) We investigate the 
quarkonium energy levels and tt threshold cross 
section in the following. 

In the phenomenological analysis we use the 
potential-subtracted (PS) mass [S] to make our 
prediction infrared renormalon free. The relation 
between the pole and PS masses is given by 

1 f d?q ~ 

m = mps{^lf) - - -^-r^V{q)\^=f,j , (11) 

where /i/ is infrared cutoff of order mag. We 
take into account non-Coulomb corrections for 
quarkonium energy levels known from literatures 
|5I6I9I1U| . So the results are complete NNNLO as 
far as energy level is concerned. 

3.1. Bottomonium masses 

Using the analytical expression for the Ci, fi, we 
obtain the relation between mass of bottomonium 
and bottom quark. It might be instructive to 
show the results using pole and PS masses: 

Mt(is) = 2mb + [1 + 1.09NLO 

-h(l.42 + 0.36„c)n.lo + (2-29 + 0.28„c)n3lo] 

= 2mb,ps + sj°^s[l + 0-19NLO 

-f (0.07 - 0.23„c)n2lo + (0-09 " ^■^^^c)^s^o] ' 
where mf, = 5 GeV, mpsi2GeV) = 4.6 GeV. 
The numbers are given separately for Coulomb 
and non-Coulomb to show numerical dominance 
of the former (in the pole scheme). One can see 
the presence of anomalously large Coulomb cor- 
rection (IR renormalon) in the pole-mass scheme, 
while in the PS mass scheme this behavior is im- 
proved and convergence of the series became bet- 
ter. 

We use the mass relation between T(15) and 
TOb_ps to extract the bottom PS mass from the 
experimental value Mx(is)\exp. = 9.460 GeV. We 
obtained 

mfc,ps(2GeV) = (4.57 ± O.OSpert. 

±0.01„, ± 0.07„o„ — pert. 

) GeV, (12) 




4.4 [ 
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Figure 1. The bottom PS mass, to&_ps(2 GeV) 
extracted from the Upsilon mass 

where the subscripts denote the source of errors. 

In Fig n we show a scale dependence of the ex- 
tracted PS mass mf,^ps(2GeV). Using extracted 
PS mass we are able to predict the masses of 
excited states of the spin triplet S-wave T fam- 
ily. In Fig|21 we show the T(2S') mass using 
mf,_p5(2GeV) — 4.57 GeV as a function of the 
renormalization scale fi. One can see that the 
large NNNLO corrections are preferable for scale 
/X > 2 GeV, however the prediction overshoots the 
experimental value My{2S) = 10.023 GeV at the 
lower scale. The naively expected natural scale 
for the bottomonium is n = Cpcusmb^ps/n (n is 
the principle quantum number), which is = 1.23 
GeV for T(2S'). This may indicate break down of 
perturbative computation for the excited T fam- 
ily. Similar behavior is observed for T(35'). 

3.2. Toponium 

In future linear colliders remnant of toponium 
15* resonance should be observed as a peak po- 
sition of the tt cross section. This enables us to 
measure the Mf^i^ig^ and extract the top quark 
mass from the data. Here we perform an exer- 
cise, how precisely we can predict the IS topo- 
nium mass when we fix the top quark mass as 
an input parameter. Adopting mt_ps = 175 GeV 
and II ~ 32.6 GeV {— Cpasmt^ps), we obtain 

Mtt{is) = (350 + 0.85 + 0.05 - 0.13 + 0.01) GeV 
= 350.78 GeV. (13) 

Scale variation between 15 < /i < 60GeV changes 
the total number only by 60 MeV. The small 
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Figure 2. T(2S') mass as a function of 
scale fj,. The lines are plotted for LO(long 
dashes), NLO (long-short dashes), NNLO(dots), 
NNNLO (solid) and shaded band shows 

Mr(2S)\exp.- 

higher order correction implies that precise top 
quark mass extraction is possible in principle, in 
the total cross section measurement. (There are 
several issues in the total cross section measure- 
ment at linear collider experiments, see for in- 
stance 22)- To obtain MS mass, being more 
commonly used in high energy processes, from 
Mtt(is) we need to know the relation between 
the PS and MS masses at 4-loop order, which is 
currently unknown. Analysis using large-Po ap- 
proximation to 4-loop MS - pole mass relation 
and direct extraction of MS mass from M^ff^igj 
is available from ref . ^2 > which is consistent with 
our result. 

3.3. The Coulomb wave function at the ori- 
gin and Green function 

In this subsection we discuss the Coulomb wave 
function and the Green function. Since the com- 
plete NNNLO non-Coulomb corrections and Wil- 
son coefficient of the production current in EFT 
are unknown, we shall discuss only the Coulomb 
corrections. As we demonstrated in the previous 
section, applicability of this method to bottomo- 
nium system is doubtful due to large NNNLO cor- 
rection and slow convergence. Thus we focus on 
the case of toponium Coulomb wave function and 
Green function. 

We show numerical formula for the IS" topo- 
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Figure 3. The Coulomb wave function at the ori- 
gin squared as a function of the scale /i. 

nium Coulomb wave function at the origin 



a.. 



-t-3.661L) 



+ a'l 



aj(5.832- 5.112L 

V 



13.73 + 6.446 In 



-22.91 + 18. ITL^)] , 



0.4333 
8.933^2) 
-) + 39.72 L 

(14) 



where L = In^fi/Cpctsfnt)- To draw Fig|3| we 
rewrote the ea. (fT^ using mt,ps{20GeV) = 175 
GeV and took into account the mass correction 
eq. (|ll|l at given order of nonrelativistic expansion 
consistently. The figure shows that perturbative 
corrections to the Coulomb wave function is rea- 
sonably small for ^ > 25 GeV and the scale de- 
pendence is mild, while in lower scale /i < 25 GeV 
the corrections are too large so the perturbative 
expansion is unreliable in the lower region. 
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Figure 4. Top quark pair production (Coulomb 
corrections only) for fi — 30 GeV, successively 
including higher order corrections. 
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In Fig^we show the threshold cross section 
as a function of CM energy for ps(20GeV) = 
175 GeV including only the Coulomb correction 
successively from LO to NNNLO. The line de- 
noted as "NNNLO exact" is a cross section ob- 
tained by numerically solving shroding equation 
for the Green function using NNNLO Coulomb 
potential. The result show a convergence of the 
perturbative approximation to the NNNLO-exact 
line. 
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Figure 5. The NNNLO tl cross section ( Coulomb 
corrections only) for ^ = 15,30,60 GeV for the 
perturbative approximation, and /i = 15 GeV for 
the NNNLO exact result. 

In FigEl we show the scale dependence of 
NNNLO cross section for the perturbative Green 
function with = 15,30,60 GeV, and /x = 15 
GeV for NNNLO-exact. The exact Green func- 
tion is stable against scale variation (so we showed 
only the case /i = 15 GeV), while perturbative 
Green function (cross section) is unstable against 
scale variation from 15 to 30 GeV (moderate 
change from = 30 to 60 GeV). This is con- 
sistent with wave function analysis, where the 
higher order corrections were large for /i < 25 
GeV, so we may conclude that the perturbative 
expansion is not reliable in lower scale region. 
The NNNLO-exact contains higher order inser- 
tions of the Coulomb potential to all order in 
eq. (|SJ . The perturbative cross section agrees well 
with NNNLO-exact for large scale where it is sup- 
posed to be reliable from the wave function anal- 
ysis. We believe that the NNNLO-exact cross 
section is reliable in wider range of /i, and the 



perturbative cross section is reliable only in the 
region /i > 25 GeV. Indeed we find that the multi- 
ple insertions of the Coulomb potential give large 
contributions to the perturbative Green function, 
and is slowly converging at small scale. Thus we 
conclude that the "correct" scale choice for the 
perturbative (Coulomb) cross section is /.t > 25, 
while choice of small scale may lead to mislead- 
ingly large uncertainties. We estimated yet un- 
known higher order Coulomb corrections should 
be less than 5 %. 
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